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The results of numerous investigators have shown the importance of labyrinth 
seal inlet swirl on the calculated dynamic stiffness of labyrinth seals. These 
results have not included any calculation of inlet leakage swirl as a function of 
geometry and sealing conditions of the given seal. This paper outlines a method of 
calculating the inlet swirl at a given seal by introducing a radial chamber to 
predict the gas swirl as it goes from the stage tip down to the seal location. 

For a centrifugal compressor, this amounts to including the flow path from the 
Impeller discharge, down the back of the disk or front of the cover, then into the 
shaft seal or eye packing, respectively. The solution includes the friction factors 
of both the disk and stationary wall with account for mass flow rate and calculation 
of radial pressure gradients by a free vortex solution. 

The results of various configurations are discussed and comparisons made to 
other published results of disk swirl. 


INTRODUCTION 

Recent reports In the literature (1,2, 3, 4, 5) have addressed the problem of 
calculating the rotordynamic coefficients for a labyrinth seal having a given inlet 
gas swirl, pressure drop, and resulting mass flow. The centrifugal design engineer 
has at his disposal from standard aerodynamic design codes the gas swirl and 
pressures at the impeller tip. The solution of the leakage gas path swirl and 
resulting pressure distribution is important not only for labyrinth seal coefficients 
but also for proper thrust balance calculations. This paper presents an approximate 
method of calculating these desired parameters using a modified version of the 
solution technique as proposed by Iwatsubo (1) and later extended by Childs and 
Scharrer (2). The extensions and modifications to the theory as outlined by (2) will 
be discussed in this paper. 

The equations of the modified formulation have been incorporated into a single 
labyrinth seal analysis computer code to permit rapid evaluation of different design 
conditions. The accuracy of the reported solution technique will be compared to 
experimental and analytical solution results as reported in reference (6) and to the 
limit case condition of zero leakage (approximately zero for computer program 
results). 


f ii£CEDir« PAGE PE**® «WT f SLMED 


225 


NOMENCLATURE 


Values are given in both SI and U.S. Customary Units. The calculations were 
made in U.S. Customary Units. 

b - tangential velocity ratio, U/ (rui) 

c* - wall separation at i tn radial chamber and radial clearance for seal tooth, m 
(in) 

- height of i ^ seal element tooth, m (in) 

- axial length of i^ seal chamber, m (in) 

m - mass flow rate, Kg/sec (lbm/sec) 

P - gas pressure, N/m (lb f /in ) 

r - radial position on disk, m (in) 

R,R q - outer radius of disk, m (in) 

R - gas constant, (Ibf * ft/lbm/°R) 

RM - mean chamber radius, m (in) 

RR - average radius of rotor seal chamber surface area, m (in) 

RS - average radius of stator seal chamber surface area, m (in) 

2 

R e - disk Reynolds number, = R^oa/v 

s - wall separation, m (in) 

S - separation ratio, = s/R Q 

2 

SJ - leakage parameter, = Vs/(R oj) 

o 

T - gas absolute temperature, °K (°R) 

U t - tangential velocity, m/s ec (in/sec) 

V r - radial velocity, m/sec (in/sec) 

Z - gas compressibility 

p - gas density, Kg/m (lbm/in ) 

t - radial wall shear stress, N/m^ (lbf/in^) 

t ^ - radial disk shear stress, N/m^ (lbf/in^) 

v - gas kinematic viscosity, m^/sec (in^/sec) 

X - position factor, = (R q - r)/R Q 

a) - rotor (disk) angular spin velocity, sec - * 

Wg - swirl velocity for free vortex, sec - * 
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METHOD OF SOLUTION 


The solution technique proposed was developed to permit the same basic theory 
and computer code to calculate both the circumferential swirl and pressure 
distribution down the disk and across the labyrinth as one coupled system. This 
technique was initially intended as a quick first pass method but the results have 
proven to be very close to the more exact theories such that the added complication 
of coupling different theories and matching boundary conditions of pressures, 
temperatures, flows, and swirls may not be justified for rotordynamic evaluations. A 
typical configuration is shown in Figure 1 for a centrifugal compressor stage disk 
cover leakage flow path. 


For the radial direction down the disk leakage path the equilibrium equation is 
given by the following equation (6): 
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dV 

r 
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For a free vortex flow neglecting the radial shear force 
by the crossflow factor in the circumferential equations 
equation becomes 

p(ru ).) 2 


or 
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which will be accounted for 
, the pressure distribution 
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where 


= gas swirl angular velocity at the radius, r. 


Therefore, the pressure at any radius r is given by 

a) 2 

P(r) = P(R) - p y-(R 2 - r 2 ) [4] 

This equation predicts the pressure along the disk if the gas swirl is known. The 
gas swirl can be calculated from the circumferential momentum equation as outlined in 
(2) and further modified to the following equation to account for variation in rotor, 
stator and mean flow chamber radius. In addition, the crossflow turbulence 
correction factors may be included in the solution to account for the inward flow 
resistance. 
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RS = mean stator surface radius 
C3, C4 = crossflow turbulence factors 
Vo^ = RM^ u)^ = average chamber swirl velocity 
DHY = hydraulic diameter of chamber 

YNR, YNS, YMR, YMS = turbulence factors per reference (2) 

ARL, ASL = shear area for rotor and stator 

p oi 

ZRT = pressure and temperature dependent gas density 
* chamber 

v = gas kinematic viscosity 

The leakage flow, m, can be calculated as outlined In (2) or by other suitable 
calculations with the modification to radial chamber pressures given by equation (4). 

The solution process requires that an Initial swirl be selected to calculate the 
pressure field and leakage. A swirl of 50% is suggested for starting the solution. 
With this pressure field and flow, the momentum equation given by Eq. (5) is solved 
for the first pass swirl values. These swirl values are then used to recalculate the 
pressure drops down the disk and through the labyrinth and the resulting flow. 
Another pass through Eq. (5) solution yields the second pass swirl values. 
Typically, three passes give the desired convergence and the inlet swirl to the laby- 
rinth is then taken from the chamber ahead of the first sealing tooth. 

A general geometry input is used such that for the radial chambers a very small 
tooth height and length with a tooth clearance equal to the disk to wall spacing can 
be used to model the flow path. The radial surface area is calculated using the 
indicated radius of each tooth location. 


in i th 


RESULTS OF SWIRL PREDICTED IN GAS LEAK PATH 


The evaluation of the proposed swirl calculation procedure has been based upon 
numerous similar conditions as reported by Jimbo (6). Initial comparisons of actual 
compressor swirl results from similar geometry is overplotted in Fig. 2. The 

parameter for leakage flow was noted to be similar to those given by Jimbo. The 

leakage flow parameter is defined as 



where 

V = radial gas velocity 
R = disk outer radius 


[ 6 ] 


a) = rotor speed 
s = wall separation 


For SJ = 0.0002 it is obvious that the compressor disk swirl does not agree with the 
reported complete analytical solution. The swirl rates are greater from the 

labyrinth program approximate solution. However, the other parameter, the disk 

Reynolds number given by 
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[ 7 ] 


2 

Re = R w/v 
o 

was calculated and found to be considerably different from the test results. The 
analysis results by Jimbo used air and a R = 9.82 x 10^ was quoted. These Initial 
compressor test results had an R = 1.1 x 10°. To match the parameters for the two 
systems the gas, pressure drop, and rotor speed were changed as given in Table 1 
under test rig conditions. This gave a disk Reynolds number of R g = 3 x 10 , only 
off by a factor of 3 from Jimbo, compared to a factor of 112 for the compressor stage 
results. 

The results for the labyrinth analysis overplot to the accuracy that the curves 
can be evaluated. The swirl results for the compressor stage and the test result 
condition are given in Tables 2 and 3, respectively. The swirl down the disk and 
through the labyrinth are given in the table with the radius /“Positl 011 factor 
indicated for comparison to the analysis results of Figure 2. The results are in 
excellent agreement for the case of near zero flow (i.e., swirl ~0.5) and for SJ = 
0.0002 where the swirl at x - P° si tion factor of 0.31 is now calculated as 0.563 as 
compared to 0.63 for the compressor gas. A comparison of the compressor swirl, test 
case calculation, Jimbo calculation, and test results reported by Jimbo are shown in 
Figure 3. The test results show a slowing of the swirl that is not predicted. The 
test rig was equiped with numerous flow and pressure measuring ports in the stator 
wall and it is very possible that the cause of the test rig result reduced swirl was 
the increased surface roughness resulting from the measuring instruments. Complex 
flow fields could also be the cause of the discrepancy and are beyond the scope of 
the present analysis. Results for reduced leakage, SJ = .000052, are given in Fig. 4 
and labyrinth analysis compressor results for SJ = .0000372 overplotted. This case 
of reduced leakage compares closely even though the disk Reynolds numbers are not 
similar. The test rig results once again show a reduced swirl ratio with great 
restriction noted in the x~POsition factor range of 0.4-0. 5. The overall trend is 
similar as concluded by Jimbo. 

One additional labyrinth program result is given in Table 4 for the condition of 
leakage from the final compressor stage through a balance piston full labyrinth. 
These results have a wall spacing that reduces as the radius reduces. A swirl of 
0.82 is predicted for this geometry and gas properties, even though the flow SJ 
parameter is 0.00021 and a uniform spacing air test result would give a swirl rate 
closer to 0.6 (see Fig. 2 for x - POsition ~0.4). 

CONCLUSIONS 

(1) The proposed approximate calculation procedure produces results that are 
acceptable for rotor dynamic evaluations of labyrinth seals. 

(2) The flow parameter, SJ, and disk Reynolds number. Re, used by Jimbo to present 
results are very useful in comparing results for different designs and give 
great insight into disk swirl behavior. 

(3) Non-uniform leak path geometry can be used to increase or decrease the swirl in 
the gas leak path. 

(4) Increased stator surface roughness will suppress the swirl due to the increased 
shear drag on the swirling flow. 
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RECOMMENDATIONS 


(1) Test evaluations using current technology flow measurement capability should be 
conducted on typical compressor and turbine disk design gas leak path 
configurations. 

(2) The importance of gas properties and actual system configurations must be 
closely evaluated. 

(3) The proposed calculation procedure can be used, with a high degree of 
confidence, for entry swirl evaluation of compressor labyrinth designs. 
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TABLE 1 SYSTEM PARAMETERS FOR ACTUAL COMPRESSOR STAGE 
DESIGN CONDITIONS AND AN ASSUMED SYSTEM TO APPROXIMATE THE 
RESULTS OF TEST CONDITIONS FROM REFERENCE 6. 


PARAMETER 

COMPRESSOR STAGE 

TEST RIG CONDITIONS 
(Ref. 6 assumed conditions) 

MW 

21.33 

25.95 

2 / 

v, m /s 

2.79 x 10 -7 

4.64 x 10 -6 

(in^/sec) 

(4.32 x 10" 4 ) 

(7.19 x 10" 3 ) 

Z 

.89 

.955 

PS, N/m 2 ) 

5.90 x 10 6 

3.43 x 10 5 

(lb/ in 2 ) 

(855) 

(49.7) 

PE, N/m 2 

5.17 x 10 6 

1.01 x 10 5 

(lb/in 2 ) 

(750) 

(14.7) 

T, °K 

331.6 

301.4 

(°R) 

(602.9) 

(548.5) 

N, Hz 

207.4 

100.0 

(RPM) 

(12566) 

(6000) 

R o» m 

.154 

.154 

(in) 

(6.05) 

(6.05) 

R e 

1.1 x 10 8 

3 x 10 6 
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TABLE 2 RESULTS FOR COMPRESSOR STAGE GIVING PREDICTED 
CHAMBER FLOW SWIRL RATIO 


CHAMBER 


SJ = 0.000023 

0.00014 

0.000236 

0.00034 

DISK TIP 

0 

.52 

.52 

.52 

.52 

2 

.08 

.5033 

.534 

.548 

.557 

3 

.17 

.5053 

.5496 

.571 

.5867 

A 

.25 

.5095 

.5689 

.600 

.6225 

5 

.31 

.5144 

.5938 

.632 

• 6603 

RADIAL t_6 

.34 

.5045 

.5997 

.6425 

.6731 

TURN 7 

- 

.452 

.5605 

.6111 

.6472 

SEAL 8 


.384 

.5121 

.571 

.6131 

+ 9 

- 

.359 

.475 

.5377 

.5834 

10 


.350 

.447 

.5097 

.5574 

11 

- 

.346 

.425 

.4862 

.5347 

12 

- 

.345 

.408 

.4663 

.5146 



R 2 oj 

x 10 8 





R = — = 1.1 
e v 



TABLE 3 

RESULTS FOR TEST CASE GIVING PREDICTED CHAMBER 


SWIRL 

RATIO 

FOR ASSUMED SUPPLY PRESSURE 

CONDITIONS 




(SEE TABLE 

1) 



CHAMBER 


SJ = .000015 

0.0001 

0.00022 

0.00034 

Disk Tip 

0 

.52 

.52 

.52 

.52 

2 

.08 

.496 

.509 

.522 

.534 

3 

.17 

.496 

.513 

.532 

.549 

4 

.25 

.497 

.519 

.544 

.567 

5 

.31 

.492 

.529 

.563 

.592 

RADIAL t 6 

.34 

.47 

.524 

.564 

.597 

TURN 7 

- 

.43 

.473 

.518 

.558 

SEAL 8 

- 

.3513 

.4105 

.465 

.511 

+ 9 

- 

.3445 

.3798 

.431 

.478 

10 

- 

.3438 

.364 

.408 

.453 

11 

- 

.3438 

.3558 

.392 

.434 

12 

- 

.3438 

.3512 

.382 

.421 


R 2 w 

o 
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TABLE 4 RESULTS OF LEAKAGE AND GAS SWIRL 
FOR FLOW FROM LAST STAGE TO A BALANCE PISTON 
LABYRINTH HAVING 15 TEETH 


CHAMBER 


X 


SWIRL 

RADIUS 

WALL SPACE 



(dim.) 


(dim. ) 

(in.) 

(in.) 

0 


0 


.637 

9 

- 

1 


.001 


.641 

8.94 

.15 

2 


.07 


.6756 

8.34 

.4 

3 


.15 


.705 

7.65 

.36 

4 


.22 


.732 

7.02 

.32 

5 


.293 


.765 

6.36 

.3 

RADIAL t 6 


.375 


.81 

5.62 

.26 

TURN 7 


.42 


.823 

5.18 

.22 

SEAL + 8 


- 


.781 

5.14 

- 

9 


- 


.74 

5.14 


10 


- 


.71 

5.14 

- 

11 


- 


.68 

5.14 

- 

12 


- 


.66 

5.14 

- 

13 


- 


.64 

5.14 

- 

14 


- 


.62 

5.14 

- 

15 


- 


.61 

5.14 

- 

16 


- 


.598 

5.14 

- 

17 


- 


.588 

5.14 

- 

18 


- 


.580 

5.14 

- 

19 


- 


.573 

v 5.14 

- 

20 


- 


.567 

5.14 


21 


- 


.563 

5.14 

- 

N 

__ 

11097. RPM 


SJ 

* 0.00021 


MW 

PS 

PE 

s= 

13.3 
948 PSI 
253 PSI 


*e 

9 2 (1162) 

• = 1.23 x 10 8 

=5 


7.63 x 10" 4 

V 

as 

7.63 x 10" 4 

in 2 / 

sec C 

= 0.52 


Z 

= 

0.979 


P 



Y 

SB 

1.255 


c 

= 0.0115 in. 

radial clearance 


leakage = 1.09 lb m /sec type seal = interlocking 
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Figure 1. - Typical compressor stage showing disk cover gas 
leakage path with nomenclature for analysis. 
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b: TANGENTIAL VELOCITY RATIO 



X- POSITION FACTOR 

Figure 2. - Velocity profile along radius, b = U/(ru) 
versus x = (R 0 -r)/R 0 at R e = Rg u/v = 9.82x105. 
Overplot of typical compressor result with 
R e - 1.1x108 (From Jimbo - fig. 6 ref. (6)). 



Figure 3. - Tangential velocity profile along radius 
b = U/(ru) versus x = R 0 -r/R o for constant 
leakage flow SJ = (s/R 0 ) [V/R 0 u) ] = 0.000104 at 
R e = 9.82x105. Overplot of a typical compressor 
result with SJ = 0.00014 and R e = l.lxlO 8 ; test 
case points for SJ = 0.0001 and R e = 3xl0 6 (over- 
plot on fig. 18 of ref. (6)). 
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